Rotating vortex dipoles in ferromagnets 
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Vort ex- anti vortex pairs are localized excitations and have been found to be spontaneously created 
in magnetic éléments. In the case that the vortex and the antivortex have opposite polarities the pair 
has a nonzero topological charge, and it behaves as a rotating vortex dipole. We find theoretically, 
and confirm numerically, the form of the energy as a function of the angular momentum of the System 
and the associated rotation frequencies. We discuss the process of annihilation of the pair which 
changes the topological charge of the System by unity while its energy is monotonically decreasing. 
Such a change in the topological charge afîects profoundly the dynamics in the magnetic System. 
We finally discuss the connection of our results with Bloch Points (BP) and the implications for BP 
dynamics. 
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I. INTRODUCTION 

Magnetic vortices have a nontrivial topology which is 
responsible for their stability and for the fact that they 
are central in theoretical studies for the micromagnetic 
description of two-dimensional (2D) magnets and thin 
films Similar topological structures arise in many 
physical Systems. In particular, magnetic vortices have 
been discussed extensively in quantum field theory 0. 
However, expérimental évidence for their existence and 
properties had been rare. This is not surprising and is 
mainly due to the fact that the vortex energy is increasing 
logarithmically with the System size and it thus diverges 
for Systems without boundaries. 

The situation has dramatically changed in the last 
years. It was realized that a disc-shaped magnetic él- 
ément, with a diameter of a few hundreds of nanometers, 
provides an excellent geometry for a magnetic vortex con- 
figuration. In particular, its exchange energy is finite 
because the System is finite, while its magnetostatic field 
vanishes everywhere except at the vortex core, thus mak- 
ing the vortex the lowest energy magnetic state. In a few 
words, the interest in the vortex stems from the fact that 
this is a nontrivial magnetic state which can, neverthe- 
less, be spontaneously created in magnetic éléments (see, 
e.g-, i). 

This leads naturally to the question whcther there are 
any further nontrivial magnetic states which would play 
an important rôle in magnetic éléments An an- 

swer comes from a somewhat unlikely direction. Récent 
experiments have shown a peculiar dynamical behaviour 
of vortices and magnetic domain walls when thèse are 
probed by external fields. Vortices may switch their po- 
larity under the influence of a very weak external mag- 
netic field of the order of a mT [6]. The same switch- 
ing phenomenon was observed by passing an a.c. elec- 
tical current in a magnetic disk [7,]. Since the polarity 
of the vortex contributes to its topological structure, the 
switching process clearly implies a discontinuons change 



of the magnetic configuration. This is certainly a surprise 
since the external field is only very weak. The key to the 
phenomenon is the appearance of vortex pairs which are 
spontaneously created in the vicinity of existing vortices 
[a, Q. The création of topological excitations (vortex 
pairs) by alternating external fields has been anticipated 
by a collective coordinate study 

We will study vortex-antivortex pairs (VA-pairs) and 
argue that thèse are nontrivial magnetic states which 
play an important rôle in dynamical phenomena in mag- 
netic cléments. They are localized configurations - unlike 
a single vortex - and they have a finite energy even when 
found in infinité Systems. Specifically, we will study a 
VA-pair where the vortex and the antivortex have oppo- 
site polarities. This behaves as a rotating vortex dipole. 
Its topology and its dynamics are radically différent than 
a pair with same polarities which was shown to undcrgo 
translational (Kelvin) motion fïd\. Despite its nontrivial 
topology a rotating vortex dipole can be destroyed by a 
quasi-continuous process. In particular, no energy barrier 
has to be overcome in contrast to the usual case for topo- 
logical solitons. This opens the possibility for switching 
mechanisms between topologically différent states in fer- 
romagnets. On account of a direct link between topology 
and dynamics in magnetic média pH IT^ , the possibility 
to change the topological magnetization structure leads 
to a dramatic change of the magnetization dynamics as 
the VA-pair is created or annihilated. Such a pair an- 
nihilation process lies at the heart of a counterintuitive 
vortex polarity switching event which was observed in 
magnetic éléments [1, 0] ■ 



II. THE MODEL 

A ferromagnet is characterized by the magnetization 
vector m — (toi, TO2, TO3) which is a function of position 
and time m = m(r, t). It has a constant length which we 
choose unity for convenience: m? = 1. The dynamics of 



the magnetization is given by the Landau-Lifshitz (LL) 
équation which can be written in a rationalizcd form 



dm 



/ = Am -Qms 63, 



(1) 



where 63 is the unit vector in the third magnétisation 
direction. The vector m is normahzed to the saturation 
magnetization Mg. Distances are measured in exchange 
length units i^x = y/AjJ^hiMY) where A is the exchange 
constant. The unit of time is tq = \I{A-k^Ms) where 7 
is the gyromagnetic ratio. Typical values are £ox ~ 5nm 
and To ~ lOps. We will use thèse values in examples 
given later in this paper. 

In the form ([T]) the LL équation has only one free 
parameter which is the quality factor Q = K/{2TrM'^), 
where K is an anisotropy constant. We will typically use 
Q = 1 in the foUowing, which corresponds to easy-plane 
anisotropy, unless stated otherwise. Eq. ([T]) is associated 
with the energy functional 



E 



Ep + En 



(2) 



where the intégration extends over an infinité plane. 
We thus consider, for simplicity, a 2D infinité Sys- 
tem. Energies are then measured in units of AnM'^l'l^. 
The anisotropy energy term in ^ models a magne- 
tocrystalline anisotropy contribution which is intrinsi- 
cally présent in materials, but it can also serve as a sim- 
plified model for the magnetostatic term in thin films. 
We defer discussion of a magnetostatic term until we de- 
rive our main results. 

A magnetic vortex is an axially symmetric solution of 
Eq. H]) of the form 

m3 = A cose(/9), mi +im2 sine(p)e^^('^-'^''), (3) 

where (p, 4>) are polar coordinates, S is an integer called 
the winding numher, A = ±1 is the vortex polarity, and 
(j)o is a constant which will be called the vortex orienta- 
tion. The magnetization angle 8(p = 0) = at the vor- 
tex center, and 0(p — > 00) = tt/2 away from the vortex 
center which corresponds to an in-plane magnetization 
at spatial infinity. The winding number S" is a topolog- 
ical invariant. The usual case for a vortex observed in 
experiments is iS = 1 and psio = i7r/2, while in the case 
S' = — 1 it is termed an antivortex. 

A further topological invariant is defined by 



A/" = / nd'^x, 
Ait ' 



^e^,y [d^m X d^m) ■ m, (4) 



where t^i, is the totally antisymmetric tensor (/i, — 
1,2), and n is a topological density. For a vortex of the 
form ^ we have N = —1/2S\. That is, for the usual 
vortices with 5 = 1 we have M — ±1/2 where the sign 
dépends on the polarity. 




FIG. 1: (colour online) VA-pairs with opposite polarities. The 
vectors show the projection of m on the (x,y) plane, ms is 
coded in colours. (a) The distance between vortices is rel- 
atively large (d = 2.7) and thus each vortex profile resem- 
bles that of an isolated one. The pair rotâtes with uj = 0.06 
[l = 68). (a) The distance between the vortices is small 
(d = 1.85) and thus the configuration is similar to Rota- 
tion frequency lj = 0.17 = 15). 



III. A VORTEX- ANTIVORTEX PAIR 

A convenient formulation for the description of vor- 
tices, but also of multi- vortex states, is obtained through 
the stereographic variable 



rrii + i 7712 

1+1713 



(5) 



A simple configuration for a two-vortex solution may be 
written as 



C 



(6) 



where a is a constant which will be considered real for 
simplicity. It gives the distance between vortices as well 
as the size of each vortex core. A vortex similar to ([3]) 
with S — l and A = — 1 is centered at C = 0,, an antivor- 
tex {S — —1) with opposite polarity A = 1 is centered at 
C = —a. At large distances |C| — > cx) we have il ^ 1. In 
other words we have chosen the boundary condition 



m 



(1,0,0) 



as 



ICI 



(7) 



which is appropriate in the presently studied model of 
an easy-plane anisotropy. The boundary condition ([7]) 
indicates that - unlike in the case of a single vortex - 
a VA-pair is a localised configuration, in the sensé that 
the magnétisation goes to the vacuum at spatial infin- 
ity. Fig.[T] shows examples of VA-pairs, whose détails are 
calculated and will be discussed later in the text. 

The form © will be used as an ansatz for a VA-pair 
with opposite polarities. It belongs to a family of exact 
static solutions of the Landau-Lifshitz équation for an 
isotropic 2D ferromagnet {Q — 0) 0, Its energy is 
E = Ee = 4:Tr for every a, (due to the scale invariance of 
the exchange interaction Q]), and its topological charge 
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is Af = 1. The topological density is even with respect 
to X n[x,y) = n{—x,y), and thus precisely half of the 
contribution to the topological charge cornes from the 
vortex in the onc half-plane {x > 0) and the second half 
cornes from the antivortex in the second half-plane {x < 
0). The vortex and the antivortex are not overlapping 
irrespectively of the distance between them. 

In the présence of an easy-plane anisotropy term {Q > 
0) the profile of the vortex and the antivortex will not 
be similar to that suggested by ansatz ^ but it will be 
modified foUowing the scale introduced by the anisotropy. 
The vortex core size is set to iîc ~ 1/ VQj ^ v^V^ 

in the usual units. Otherwise, the comments on the topo- 
logical density and charge of the VA-pair given in the 
previous paragraph remain valid. Finally, the boundary 
condition ([7]) is consistent with easy-plane anisotropy. 

IV. VORTEX-ANTIVORTEX DYNAMICS 

A more detailed study of the VA-pair requires infor- 
mation on its dynamics. For this purpose we write the 
linear momentum associated with the LL équation ^ 



p. = 



y n dxdy, 



Py = 



X n dxdy, 



(8) 



whose components are actually the moments of the topo- 
logical density n. Since the vortex and the antivortex 
have a mirror image topological density distribution we 
find that {Px,Py) — (0,0). Formulae ^ give a measure 
of the position of the VA-pair (which is at the origin). 
The important resuit is that this mean position is con- 
served in time since (HJ are naturally conserved quan- 
tities. We thus see that the VA-pair is spontaneously 
pinned at the position where it is created. This dynam- 
ics should be contrasted to the dynamics of a VA-pair 
with same polarities. The latter forms a solitary wave 
which is necessarily non-static and undergoes transla- 
tional (Kelvin) motion ^ï^. 

Further information is obtained by considering the an- 
gular momentum of the System [Tî| 



p ndxdy, 



(9) 



which is also written with the aid of the topological den- 
sity n, and it is conserved within model ([T]). For a VA- 
pair £ is clearly nonzero (the integrand is positive def- 
inite), and it gives a measure of its size. The distance 
d between the vortex and the antivortex can be defined 
from 



J n dxdy i 



J n dxdy 2t: M 



(10) 



where we assume that the VA is centered at the origin. 
Using the value TV = 1 for the vortex pair we find 

d^ = -e. (11) 



A nonzero angular momentum £ indicates that the VA 
should be a rotating object. 

The rotational dynamics of the VA-pair is fuUy con- 
firmed by numerical simulations of the LL Eq. ([T]) . How- 
ever, the rotating pair is apparently unstable and short 
lived due to radiation of energy in analogy to a rotating 
electric dipole. The conservative System will be further 
studied here in order to obtain some useful results. Thèse 
will then be used to understand the full dynamics of a 
more realistic System including mechanisms for energy 
dissipation. 

For a simpler numerical investigation of the main fea- 
tures of a rotating pair, we consider a steady state of a 
rotating VA-pair for every value of the conserved angular 
momentum i. One should note that this is a stationary 
point of the extended energy functional 



F = E-uj£, 



(12) 



where w is the angular frequency of rotation of the VA- 
pair measured in units of A-k^Ms. Using the standard 
scaling argument [3] that a stationary state is a mini- 
mum of the functional with respect to scaling the 
coordinates, one readily finds, for a stationary state, the 
virial relation 



(13) 



The exchange energy does not appear in the above rela- 
tion due to its scale invariance in two dimensions. 

A useful resuit is now obtained if we assume well sepa- 
rated vortex and antivortex where each of them behaves 
as an isolated one. The anisotropy energy for an isolated 
single static vortex in the model ^ is n/2 [ï^. Substi- 
tuting this number in Eq. (|13p we obtain 



1 



2 

d2' 



(14) 



where Eq. (|lip was used in the second equality. The 
limit of large séparation between vortices was studied in 
Réf. Q through a collective coordinate approach and the 
angular frequency (fT4|) was obtained. Eq. (fT4|) shows that 
the rotation frequency goes to zéro in the limit of large 
vortex séparation. As a concrète example, suppose that 
the vortices are a distance d = lO^ex apart and Q = l, 
to obtain uj = 0.01, or a; ~ 10^ Hz. 

We now turn to the case that the vortex and antivortex 
are close and they are forming a small pair. When the 
size of the VA-pair is smaller than the length scale 1 / y/Q 
introduced by the anisotropy energy, the latter is negligi- 
ble. Then ^ is exact in the area around the origin close 
to the vortex cores. Away from the vortex cores, we will 
adopt the approximation that the anisotropy will prevail 
and it will enforce the vacuum m = (1, 0, 0). In the limit 
a ^ 1 we have E — Ee — in [ï^ and A/" = 1 due to the 
contribution of the origin which becomes a singular point. 
The magnetic configuration becomes then fuUy aligncd to 
the x-axis, except for the origin as a — > 0. Then, one can 
eliminate the singular point and assume a fuUy aligned 
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FIG. 2: The energy i5 of a rotating VA-pair as a function of 
its angular momentum i (for Q = 1). The circles indicate 
numerical results (the dashed line is a guide to the eye). It 
can be fitted by _E = 7rln(^/4.5) + 47r for large £, while E — 
0.5^ + 47r for — + 0. The rotation frequency is the slope of 
the curve. The energy of a vanishingly small VA-pair is 47r. 



State. In doing so, one éliminâtes the topological com- 
plexity of the System and changes the topological charge 
M by unity. Needless to say, such singular points may not 
exist in magnetic materials because the lattice spacing in 
the solid provides a natural length cut-ofï. The finite 
lattice spacing is apparently responsible for the annihi- 
lation of vanishingly small VA-pairs and the subséquent 
surprising change of the topological charge of magnetic 
configurations in récent experiments 0, 0| • 

One usually expects that a topologically nontrivial 
state {M ^ 0) cannot be deformed continuously to the 
ground state of the System which has topological charge 
zéro. In addition, the two states are separated by an in- 
finitely high energy barrier. However, in the case of the 
VA-pair studied here, although the former argument is 
correct, no energy barrier is encountered in the process. 
This is an unusual property but it is completely explained 
by the scale invariance of the exchange energy 0. 

We go on to find numerically VA-pairs in steady rota- 
tion. We set up a numerical grid which goes to spatial 
infinity by using stretched coordinates ry = a tan a;, ^ = 
6 tan y. The choice of the constants a, 6 dépends on the 
size of the studied vortex pair and they are chosen so 
that it is described at sufîicient resolution. We first go to 
a rotating référence frame which amounts to substitut- 
ing /—*■/ — u){5t/5m) in Eq. We use a relaxation 
algorithm, i.e., we add dissipation to the LL équation 
(effectively using Eq. p7|) ). We feed the algorithm with 
the ansatz ([S]) as an initial condition, and this subse- 
quently relaxes to a roughly steady state. Our algorithm 
cannot give exact solutions of the équation; a more elab- 



orate numerical method would have to be used for that 
purpose. However, using appropriate parameters for the 
numerical lattice, we are able to find approximate solu- 
tions which are almost steadily rotating states. We thus 
find good approximations for the energy and the angu- 
lar momentum of the rotating state for varions values of 
the rotation frequency lo. Fig.[T] shows rotating VA-pair 
solutions. Fig. [5] shows the results for the energy as a 
function of the angular momentum for Q = 1. It is clear 
that the energy of a vanishingly small VA-pair (at € — > 0) 
is pure exchange energy and it is equal to 47r, as antic- 
ipated. We have checked that the virial relation (fTH]) is 
approximately satisfied for our numerical solutions. We 
note that it appeared more difficult to obtain reasonable 
convergence of the algorithm for smaller values of (.. 

We can find the rotation frequency of a vanishingly 
small VA-pair by substituting Eq. ^ in Eqs. ^ and 
to find 



{27rpdp), (15) 



where S dénotes a circular domain which includes the 
vortex cores. Substituting ([TÏÏ|) in the virial relation (fT5|) 
we find lû — Q/2. A more rigorous calculation would 
require an asymptotic analysis of the LL équation. On 
the side of small £ in Fig. [21 we have a steep increase of 
u! consistent with the theoretical value u!{£ ^ 0) = 1/2 
at (Q = 1). 

We can give an analytic estimation of the curve in 
Fig. [2] for large values of the angular momentum £. We 
first note that the frequency satisfies lu — dE/d£ and it 
is thus given by the slope of the curve. For large £ we 
can substitute Eq. for bj{£) and obtain a differential 
équation which givcs 



E^TT ln(£/4) + 47r, 



(16) 



where the constant £q cannot be fixed by the présent 
calculation (the constant 47r has been added to facilitate 
comparison with Fig. [5]). The logarithmic behaviour of 
the energy can be obtained by considering the energy 
of two widely separated vortices. The anisotropy energy 
should be independent of £ while the interaction between 
the vortices is due to exchange. The exchange energy 
is approximately equal to 27r \Ti{d/ Rc) and this coïncides 
with Eq. ([Tïï]) when we use Eq. (fH)) for large £. 

We now turn to discuss more realistic Systems where 
dissipation is présent. The LL équation including Gilbert 
damping reads 

dm 

— = m X f - am X {mx f), (17) 

where a is the dissipation constant. The dissipation is 
necessary in order that the simulated results are doser to 
the expérimental, but it also helps to rid the numerical 
results of any fast transients and thus focus on the main 
features. In particular, we typically consider in this pa- 
per strong enough Gilbert damping so that it dominâtes 



5 



the radiation efFect disussed earlier in connection to con- 
servative dynamics. 

We typically use configuration ([6]) as an initial state for 
the numerical simulation. This is iterated in time using, 
typically, a dissipation constant a ~ 0.1 and values of a ~ 
5. In the initial phase of the simulation the vortex and the 
antivortex cores quickly shrink to the size Rc discussed 
above. Then the vortex and the antivortex approach each 
other and thus the size of the pair decreases. This is 
certainly a resuit of dissipative dynamics (the Gilbert 
term in Eq. (|17p ) as the energy (mainly the exchange 
part) is obviously decreasing with the vortex pair size. 
The rotational dynamics (due to the conservative term 
in Eq. (fT7|) ) is at the same time présent and it becomes 
prévalent for small dissipation a. 

The dynamical behaviour of a VA-pair can now be 
summarized as follows. A vortex and an antivortex ro- 
tate around each other, while the pair shrinks due to 
dissipation of energy. Its energy follows approximately 
the curve of Fig. [2] as its size (and its angular momen- 
tum) decreases. At vanishing size a singular point, as 
the one discussed earlier, would be created. This would 
practically disappear due to the discreteness of the solid. 

We now turn to discuss the relevance of the above ideas 
in a thin film environment. The addition of the magneto- 
static field in Eq. (P) or Eq. ([TT]) is now necessary in order 
to make the model relevant for experiments in magnets. 
The qualitative picture for the dynamics of a VA-pair 
will not be affected since our analysis does not dépend 
on the exact form of the field / in Eq. ([T]). The mag- 
netostatic field will, however, modify the profile of the 
pair. Unlike the exchange and anisotropy, the magneto- 
static energy is sensitive to rotations of the spin vector. 
It has been pointed out @ that this symmetry breaking 
is responsible for the résonant excitation of vortex pair 
by alternating external fields. 

For an ultra-thin film, we can assume that there is no 
variation of m across the film thickness, and we then 
expect that the magnetostatic interaction would behave 
as a single-ion easy-plane anisotropy (as for the ground 
state in thin films [13] )■ This amounts to the substitu- 
tion Q ^ Q + 1 in Eq. ([T]). Our preceeding analysis 
is then valid for an effective anisotropy strength Q + l. 
In particular, a change in the topological charge occurs 
through a 2D singular point. It is interesting to note that 
in this particular case no Bloch Point need be formed. In 
the spécial case of a significant easy-axis magnetocrys- 
talline anisotropy, at Q = — 1, we effectively recover the 
isotropic model. The curve in Fig. [2] would then be a 
straight line E{£) = An. Consequently, a VA-pair which 
is somehow created in the magnet would be static. More 
realistically (experimentally) one would expect that it 
would be slowly moving. 

It would be interesting to consider the behaviour of 
the curve as the thickness of the film is increased. A 
fuU study can only be donc numerically and will not be 
given here. An interesting possibility is that the long- 
range magnetostatic interaction would compete with the 



local anisotropy interaction to give a non-monotonie E{i) 
curve, for values of Q ~ — 1- A minimum in the curve 
for small £ could then appear because the vortex and the 
antivortex have opposite polarity which may resuit in a 
low magnetostatic energy. If such a minimum is there, 
we would have a; = at that point and thus an unusual 
state of a static VA-pair would be stabilized. 

Numerical simulations [1, [ï^ show that, in a film of 
finite thickness, a singular point is created first at one 
of the surfaces of the film. This should apparently be 
attributed to the variation of the magnetization across 
the film thickness. The VA-pair vanishes by formation 
and annihilation of the singular point at successive film 
levels, until the vacuum state is reached throughout the 
film. A curve similar to that in Fig. [5] gives now the 
linear energy density across the film thickness. The en- 
ergy released when the VA-pair is annihilated is Ant (t is 
the film thickness) and it is emitted through spin waves 

At the stage when the singular point has been formed 
and annihilated, say, near the top film surface, while a 
VA-pair still exists in the lower film levels, we have a 
Bloch Point (BP) in the film [Î9[. This lies at the inter- 
face between the région where the VA has been annihi- 
lated and the région where this has a finite size, at the 
singular point. This is a somewhat simplified realization 
of the BP studied in Réf. [1^. Its energy is approxi- 
mately An times the length of the vortex or the antivor- 
tex line. More important is the unusual fact that during 
création and annihilation of the BP one does not have to 
overcome any energy barrier (unlike the case discussed 
in Réf. [13]), essentially for the same reasons explained 
above in connection to the singular point. 



V. CONCLUSIONS 

We have studied a vortex antivortex pair where the 
two vortices have opposite polarity. The pair has a 
nonzero topological number A/" = 1, and the vortices ro- 
tate around each other in analogy to a rotating electric 
dipole. As energy is dissipated the pair shrinks and it 
eventually dégénérâtes to a singular point. Such a sin- 
gular point would disappear in a solid due to the dis- 
creteness of the lattice. The latter process changes the 
topological number of the magnetic configuration. It is 
surprising that no energy barrier needs to be overcome 
during this quasi-continuous process. 

The mechanism for changing the topological number 
of the magnetization configuration makes it possible to 
obtain controUed switching between topologically dis- 
tinct (and therefore robust) magnetic states. This was 
achieved in the experiments of Refs. [1, The com- 
plète mechanism, which is common in both experiments 
as regards vortex dynamics, includes a dynamical pro- 
cess which involves three vortices @. Our analysis ap- 
plies only to the final part of this process, that is, when 
the vortex and antivortex with opposite polarities are 
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relatively close together. This is one of the important 
stages of the process at which the topology of the mag- 
netic configuration actually changes. The présence of a 
third vortex does not modify our main arguments. The 
fuU three-vortex System has a nonvanishing topological 
charge and thus it is a rotating object which is sponta- 
neously pinned in the magnet. A study of this configura- 
tion could be donc foUowing the methods of this paper. 

Our discussion clearly suggests a connection between 
the VA-pair dynamics and the dynamics of a BP. In par- 
ticular, a rotating motion of the magnetization appears 
inhérent in a BP configuration. This effect is apparently 



not présent in the spécial spherically symmetric BP so- 
lution [20|]. It is not présent either in the dissipative 
dynamics of a BP discussed in the literature [2^ . 
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